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ABSTRACT 

Surveys above lOkeV represent one of the the best resources to provide an 
unbiased census of the population of Active Galactic Nuclei (AGN). We present 
the results of 60 months of observation of the hard X-ray sky with Swift/BAT. 
In this timeframe, BAT detected (in the 15-55 keV band) 720 sources in an all- 
sky survey of which 428 are associated with AGN, most of which are nearby. 
Our sample has negligible incompleteness and statistics a factor of ~2 larger 
'"q_j over similarly complete sets of AGN. Our sample contains (at least) 15 bona-fide 

Compton-thick AGN and 3 likely candidates. Compton-thick AGN represent a 
~5 % of AGN samples detected above 15 keV. We use the BAT dataset to refine 
the determination of the LogN-LogS of AGN which is extremely important, now 
that NuSTAR prepares for launch, towards assessing the AGN contribution to the 
cosmic X-ray background. We show that the LogN-LogS of AGN selected above 
10 keV is now established to a ~10 % precision. We derive the luminosity function 
of Compton-thick AGN and measure a space density of 7.9l|^ x 10"^ Mpc'^ for 
objects with a de-absorbed luminosity larger than 2x10^^ erg s"^. As the BAT 
AGN are all mostly local, they allow us to investigate the spatial distribution of 
AGN in the nearby Universe regardless of absorption. We find concentrations 
of AGN that coincide spatially with the largest congregations of matter in the 
local (<85Mpc) Universe. There is some evidence that the fraction of Seyfert 2 
objects is larger than average in the direction of these dense regions. 
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X-rays: diffuse background - surveys 
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Introduction 



There is a gen eral consensus that t he cosmic X-ray background (CXB), discovered more 
than 40 years ago (iGiacconi et aLlll962l ). is produced by integrated emission of Active Galac- 
tic Nuclei (AGN). Indeed, below ~3keV sensitive observations with Chandra and XMM- 



Newton have d irectly resolved as much as 80% of the CXB into AGN (jWorsley et al.l 12005 



Luo et al.ll201ll ). However, above 5 keV, due to the lack of sensitive observations, most of the 
CXB emission is at present unresolved. Po pulation synthesi s models have successfully shown, 
in the context of the AGN unified theory (lAntonuccil 119931 ). that AGN with various level of 
obscuration and at dif f erent redshifts can account for 80-100% of the CXB up to ~100keV 
( IComastri et al.l Il995t iGilli et al.l l2001t iTreister fc Urryl 120051 ). In order to reproduce the 
spectral shape and the intensity of the CXB, these models require that Compton-thick AGN 
(Nh > 1.4 X lO^^cm"^) contribute ~10% of the total CXB intensity. With such heavy 
absorption Compton-thick AGN have necessarily to be nur nerous, comprising p erhaps up to 



30-50 % of the AGN population in the local Universe (e.g. iRisaliti et al.lll999[ ). However, it 
is still surprising that only a very small fraction of the populatioii of Compton-thick AGN 
has been uncovered so far (IComastrill2004j : iDella Ceca et al.ll2008al and references therein). 



Studies of AGN are best done above 10 keV where the nuclear radiation pierces through 
the torus for all but the largest colur nn densities. Focusing optics like those rn ounted on 
NuSTAR and ASTRO-H (respectively, iHarrison et al.ll2010l : iTakahashi et al.ll2010l ) will allow 
us to reach, for the first time, sensitivities <10~^^erg cm~^ s~^ above lOkeV permitting us 
to resolve a substantial fraction of the CXB emission in this band. Given their small field 
of views (FOVs) those instruments will need large exposures in order to gather reasonably 
large AGN samples. Because of their good sensitivity the AGN detected by NuSTAR and 
ASTRO-H should be at redshift ~1, but, due to the small area surveyed, very few if any will 
be at much lower redshift. 

All-sky surveys, like those performed by Swift/BAT and INTEGRAL above lOkeV are 
very effective in making a census of nearby AGN, thus providing a natural extension to more 
sensitive (but with a narrower FOV) missions. Here we report on the all-sky sample of AGN 
detected by BAT in 60 months of exposure. Our sample comprises 428 AGN detected in the 
whole sky and represents a fact or of ^2 improveme nt in number statistics when compared to 
previous complete samples (e.g. lBurlon et al.ll201l[ ). In this paper, we present the sample and 
refine the determination of the source count distribution and of the luminosity function of 
AGN. This is especially important considering the upcoming launch of NuSTAR (scheduled 
for March 2012) as it allows us to make accurate predictions for the expected space densities 
of distant AGN. We also use the BAT sample to investigate the spatial distribution of AGN 
in the local Universe. We leave for an upcoming publication the follow-up of all new sources 
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using 2-10 keV data and the determination of the absorption distribution. 

This paper is organized as follows: the BAT observations are discussed in § |21 while § 13.11 
and § 13.21 discuss respectively the source count distribution and the luminosity function of 
AGN. In §|l]we present a measurement of the over-density of AGN in the local Universe, while 
in § [5] the prospects for the detection of AGN by NuSTAR are discussed in the framework 
of the BAT observations and population synthesis models. Finally, § |6] summarizes our 
findings. Throughout this paper, we assume a standard concordance cosmology (Ho=71 km 
s"i Mpc-i, fiA/=l-^]A=0.27). 



2. Properties of the Sample 



The Burst Aler t Telescope (BAT; iBarthelmy et al.l |2005| ) onboard the Swift satellite 
( iGehrels et al.l 120041 ). represents a major improvement in sensitivity for imaging of the hard 
X-ray sky. BAT is a coded mask telescope with a wide field of view (FOV, 120° x 90° par- 
tially coded) aperture sensitive in the 15-200 keV range. Thanks to its wide FOV and 
its pointing strategy, BAT monitors continuously up to 80% of the sky every day achiev- 
i ng, after several years, deep exposures across the entire sky. Results of the BAT survey 
dMarkwardt et all boos l: lAiello et al.l boOSat iTueller et all boosl ) show that BAT reaches a 
sensitivity of ~1 mCraqj in 1 Ms of exposure. Given its sensitivity and the large exposure al- 
ready accumulated in the whole sky, BAT is an excellent instrument for studying populations 
whose emission is faint in hard X-rays. 

For the analysis presented here we use 60 months of Swift/BAT observations taken be- 
tween March 2005 and Mar ch 2010. Dat a screen ing a nd processing was pe rformed according 
to the recipes presented in lAjello et al.l (l2008al ) and lAjello et al.l (l2008bf ). The chosen en- 
ergy interval is 15-55 keV. The all-sky image is obtained as the weighted average of all the 
shorter observations. The final image shows a Gaussian normal noise and we identified 
source candidates as those excesses with a signal-to-noise (S/N) ratio > 5 a. The final sam- 
ple comprises 720 sources detected all-sky. Identifi cation of thes e obje ct s was performed 



by cro s s-correlating our cata log wi th the catalog s of iTueller et al.l ( 12 0081 ). ICusumano et al. 



(120101 ) ■ IVoss fc Aiellol fl2010h . and iBurlon et al.l fcoill). When e ver available we use d the 
newest optica l ident ifications provided by iMasetti et al.l ( 120081 ) . iMasetti et al.l (120091 ) . and 
Masetti et al.l (120101 ). Of the 720 all-sky sources only 37 (i.e. ~5 %) do not have a firm identi- 
fication. This small incompleteness does not change when excluding or including the Galactic 
plane. Of the 720 objects, 428 are identified with AGN. This represents an improvement 



^ImCrab in the 15-55 keV band corresponds to 1.27x10 erg cm ^ s 
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of a factor >2 in the number of detected A G N with respect to previ ous complete samples 



Ajello et al. 


2009b: 


Burlon et al. 


2011 


)■ 



[e.g. lAjello et al.ll2009bl : IBurlon et al.l 1201 11 1. iCusumano et al.l ( l2010l ) recently reported on 
the sample of sources detected by BAT in 58 months of observations. Their catalog is con- 
structed using three energy bands and selecting >4.8(j excesses in any of the three bands. 
As such their catalog is larger than the one presented here. Howe ver, for the scope o f this 
and future analyses (e.g. a follow-up work of that presented in iBurlon et al.ll201ll ) it is 
important to have a clean sample whose selection effects are well understood and can be 
accounted for during the analysis. 

Fig. [1] shows the sky coverage of the BAT survey. It is apparent that the limiting flux 
is ~0.45mCrab (~5.5xl0~^^ erg cm~^ s~^) and that the BAT survey becomes complete (for 
the whole sky) for source fluxes >lmCrab. The sensitivity scales nicely with the inverse 
of the square root of the exposure tir ne as testified by th e limiting sensitivity of 0.6mCrab 
reached in 36 months of observations (lAjello et al.ll2009al ). 
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Fig. 1. — Sky coverage of the BAT survey for the 15-55 keV band and for sources detected 
all-sky above the 5 a level. 



2.1. Jet-dominated and Disk-dominated Objects 



Jet-dominated AG N (radio galaxies a nd blazars) constitute a ~15 % fraction of the 
BAT samples (e.g. see lAjello et al.l l2009bl ). This is confirmed also here where 67 (out of 



Table 1. The 428 AGN detected by BAT ^ 



SWIFT NAME 


R.A. 


Dccl. 


Position Error 


Flux 


S/N 


ID 


Type^ 


Redshift 


Photon Index 


Log Lx 




(J2000) 


(J2000) 


(arcmin) 


(10-" cgs) 














J0004.2+7018 


1.050 


70.300 


6.551 


0.76 


5.5 


2MASX ,100040192+7019185 


AGN 


0.0960 


2.04+0.46 


44.2 


J0006.2+2010 


1.571 


20.168 


3.276 


1.06 


6.5 


Mrk 335 


Syl 


0.0254 


2.60+0.32 


43.2 


JOOlO.4+1056 


2.622 


10.947 


2.281 


1.85 


11.0 


QSO B0007+107 


BLAZAR 


0.0893 


2.23+0.20 


44.6 


J0018.9+8135 


4.732 


81.592 


4.720 


0.94 


6.5 


QSO J0017+8135 


BLAZAR 


3.3600 


2.51+0.52 


48.3 


J0021.2-1908 


5.300 


-19.150 


4.773 


0.92 


5.1 


1RXSJ002108.1-190950 


AGN 


0.0950 


1.96+0.45 


44.3 


J0025.0+6826 


6.264 


68.436 


4.721 


0.78 


5.6 


IGR J00256+6821 


Sy2 


0.0120 


1.66+0.33 


42.4 


J0033. 4+6125 


8.351 


61.431 


4.274 


1.01 


7.3 


ICR J00335+6126 


AGN 


0.1050 


2.46+0.28 


44.5 


J0034.6-0423 


8.651 


-4.400 


6.165 


0.92 


5.2 


2MASX J00343284-0424117 


AGN 


0.0000 


1.79+0.43 




J0035. 8+5951 


8.965 


59.852 


2.095 


2.05 


14.8 


lES 0033+59.5 


BLAZAR 


0.0860 


2.74+0.18 


44.6 


J0038.5+2336 


9.648 


23.600 


5.132 


1.01 


6.2 


Mrk 344 


AGN 


0.0240 


1.80+0.59 


43.1 


J0042.8-2332 


10.701 


-23.548 


3.068 


2.52 


14.7 


NGC 235A 


Sy2 


0.0222 


1.90+0.11 


43.4 



''The full tabic is available in the online version of the paper. 
''AGN are sources lacking an exact optical classification. 
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Fig. 2. — Position on the luminosity-redshift plane of the 428 AGN detected by BAT in the 
15-55 keV band. The color coding reflects the optical classification reported in Tab. [U AGN 
are sources lacking an exact optical classification. The black squares mark the position of 
the Compton-thick AGN reported in Tab. [2j Note that their luminosities were not corrected 
for absorption (see text for details). The dashed line shows the flux limit of the BAT survey 
of 5.5x10^^^ erg cm"^ s"^. 

the 428 AGN) are classified as either radio-galaxies or blazars. The remaining 361 AGN are 
associated with objects optically classified as Seyfert galaxies (323 objects) or with nearby 
galaxies (38 sources), through the detection of a soft X-ray counterpart, for which an optical 
classification is not yet available. The full sample is reported in Tab. [TJ For all the sources, 
/c— corrected Lx luminosities were computed according to: 

where Fx is the X-ray energy flux in the 15-55 keV band and Tx is the photon index. This 
assumes that the source spectra are adequately well described by a power law in the 15- 
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55k eV band in a g reeme nt with what found by lAjello et aLl (j2008bl ). iTueller et aLl ( 120081 ). 
and iBurlon et aLl ( 1201 ll ). Unless noted otherwise (i.e. § 13. 3p luminosities are not corrected 
for absorption along the line of sight since this correction is different than uni ty (in the 
15-55 keV band) only for Compton-thick AGN (see Fig. 11 in lBurlon et al.ll201ll ) and does 
not introduce any apparent bias in any of the results shown in the next sections. 

Fig. [2] shows the position of the 428 sources in the luminosity-redshift plane for the 
different optical classifications reported in Tab. [TJ The BAT AGN sample spans almost 8 
decades in luminosity and includes sources detected from z^^O.OOl (i.e. ~4Mpc) up to z^4. 
It is also evident that Seyfert galaxies dominate the low-luminosity part of the sample, while 
blazars and radio-galaxies dominate the high-luminosity part of the sample. The increased 
exposure of BAT allows us to detect fainter AGN with respect to previous samples. Indeed, 
the average flux of th e Seyfert-like AGN d ecreased by ~20 % when comparing it to the sample 
of AGN reported in iBurlon et al.l ( 1201 ll ). Since the average redshift in the two samples is 
very similar, this translates into a larger number of low-luminosity AGN. 



2.2. Compton-thick AGN 



Hard X-ray selected samples are among the best resources to uncover Compton-thick 
AGN which are otherwise difficult to detect. A detailed measurement of the absorbing col- 
umn density of all the AGN in this sample is beyond the scope of this paper and left for a fu- 
ture publication. However, in order to determine the likely candidates, it is possible to cross- 
correlate our source list with catalogs o f Compton-th i ck AG N. Our AGN catalog contains all 
the 9 Compton-thick AGN reported bv lBurlon et al.l (120111) and 6 a d ditiona l Compton-thick 



AGN reported in the list of bona-fide objects of iDella Ceca et al.l (l2008a[). There 3 addi- 
tional sources which are labeled as Compton-thick candidates by iDella Ceca et al.l ( l2008al ) 
(see their Table 2) which are also detected in this sample. The full list of 18 known Compton- 
thick AGN contained in this sample is reported in Tab. [2l It is clear that the number o f 
(likely) Compton-thick AGN is doubled with respect to the sample of iBurlon et al.l (l201l[ ) 
and that Compton-thick AGN represent a 'steady' 5% fraction (i.e. ~18/361) of AGN 
samples selected above lOkeV. 

The redshift distribution of Compton-thick AGN is also different than that of the whole 
AGN sample. The median redshift of the Compton-thick AGN of Tab. |2]is 0.010 while that 
one of the entire AGN sample is 0.029. Compton-thick AGN can be detected by BAT only 
within a distance of ~100Mpc beyond which the strong flux suppression caused by the 
Compton-thick medium limits the capability of BAT to detect these objects. 
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We also checked if any of th e remaining 3 bona-fide Compton-thick AGN (or the 20 
remaining candidates) reported in iDella Ceca et aLl (j2008a[ ) he just below the reliable BAT 
detection threshold. None of the remaining sources in the above lists exhibits a significance 
larger than 3.5 a in our analysis. This means that none of these sources are likely to be 
detectable by BAT in a deeper survey. The main consequence is however that the new 
Compton-thick objects that will appear in the BAT samples will be new (i.e. previously 
un-studied) sources. A few might already be present in this sample and this aspect will be 
investigated in a follow-up study. 



- 9 - 



Table 2. Known Compton-thick AGN detected in the BAT sa mple . Unless \ yritte n 
explicitly the values of the absorbing column density come from lBurlon et al.l (120 111 ). 



NAME 


Type Reshift 


R.A. 


Decl. 








(J2000) 


(J2000) 


(1024 cm-2) 



NGC 424 


Sy2 





011588 


17.8799 


-38.0944 


1.99 


NGC 1068 


Sy2 





003787 


40.7580 


-0.0095 


>10 


NGC 1365^-'' 


Syl.8 





005460 


53.4442 


-36.1292 


3.98 


CGCG 420-015 


Sy2 





029621 


73.3804 


4.0600 


1.46 


SWIFT J0601. 9-8636 


Sy2 





006384 


91.1972 


-86.6245 


1.01 


Mrk 3 


Sy2 





013509 


93.9722 


71.0311 


1.27^ 


UGC 4203^'^ 


Sy2 





013501 


121.0585 


5.1217 


>1.00^ 


NGC 3079 


Sy2 





003720 


150.4701 


55.6978 


5.40 


NGC 3281 


Sy2 





010674 


157.9743 


-34.8571 


1.96^ 


NGC 3393 


Sy2 





012500 


162.1000 


-25.1539 


4.50 


NGC 4939 


Syl 





010374 


196.1000 


-10.3000 


>10<^ 


NGC 4945 


Sy2 





001878 


196.3726 


-49.4742 


2.20*^ 


Circinus Galaxy 


Sy2 





001447 


213.3828 


-65.3389 


4.30*^ 


NGC 5728 


Sy2 





009467 


220.6916 


-17.2326 


1.0 


ESQ 138-1 


Sy2 





009182 


253.0085 


-59.2386 


1.5^'-^ 


NGC 6240 


Sy2 





024480 


253.3481 


2.3999 


1.83 


NGC 6552^-'^ 


Sy2 





026550 


270.0981 


66.6000 


>1.00^ 


NGC 7582 


Sy2 





005253 


349.6106 


-42.3512 


1.10 



'^Par t of the sample of candidate Compton-thick objects in iDella Ceca et al. 
J2008ah . 

'^NGC 1365 is a complex source that shows a col umn density that can vary 
from Log Nff ~ 23 to > 24 on t imescales of ~10hr (IRisaliti et al.l l2009bl ) . Ac- 
cording to lRisaliti et al.l ( l2009a( ) the source has an absorber with LogNj^ ^ 24.6 
which covers ~80 % of the source. 

■^UGC 4203 (also called the 'Phoenix' galaxy) is known to exhibit changes 
in the absorbin g column density fr om the Compton-thin to the Compton-thick 
regime (see e.g. IRisaliti et al.ll2010l ). 



^Reported to be Compton-thick by lReynolds et al.l (119941 ). and lBassani et al. 



(119991). 
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Tor the value of the absorbing column density see iDella Ceca et al.l (l2008aj) 
and references therein. 



Piconcelli et al.l ( 120111 ) reports that this source might be absorbed by 
LogNn >25. 
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3. Statistical Properties 



3.1. The Source Count Distribution 



The source count distribution of radio-quiet A GN (also called LogN-LogS) has already 
been derived above 10 keV by several au thors (e.g. 



Ajello et al. 


2008b: 


Tueller et al. 


2008 



Ajello et al.ll2009bt iKrivonos et al.ll2010[ ). Here we use our larger complete set of AGN to 
refine the determination of the LogN-LogS which is of particular interest since NuSTAR will 
be surveying, with a factor >100 better sensitivity, a similar energy band in the very near 
future. Our aim is to also to compare the BAT observations to the predictions of popular 
population synthesis models. 

In order to account robustly for uncertainties in the determination of the LogN-LogS 
we perform a bootstrap analysis, creating 1000 resampled set extrac ted (with replacern ent) 
from the BAT dataset. We perform a maximum likelihood fit (see lAjello et al.ll2009bl . for 
details) to each data set with a power law of the form: 



§ = ^ * (S/IO- 



11 



(2) 



where S is the source flux, and a and A are respectively the slope and the normalization of 
the power law. From the 1000 realizations of the BAT AGN set we derive the distributions 
of the normalization and of the slope and we use these to determine the best-fit parameters 
and their associated errors. 

From our analysis we find the best-fit values of: a=2.49lo;o7 and A=1.05to;o4 x 10^. So 
the BAT LogN-LogS is compatible with Euclidean for all fluxes spanned by this analysis as 
shown in Fig. |3l The surface density of AGN at fluxes (15-55 keV) grater than 10~^^ erg cm~^ 



-1 is 6.671 ^-1^x10- 



'de£ 



reported in I Ajello et al.l (l2009al ) . 



^ whic h is in good agreement with the value of 6.7±0.4 x 10 deg 



In order to compare our results with the LogN-LogS measurements published elsewhere 
we adopt the following two strategies to convert fluxes from one band to another. In the 
flrst case we adopt a simple power law with a photon inde x of 2.0 which is kno wn to describe 
generally well the spectra of faint AGN in the BAT band (I Ajello et al.ll2008bl ). Additionally 
we use a more complex model for the AGN emission in the BAT band whic h is based on 
the PEXRAV model of iMagdziarz fc Zdziarskil Jl995h . In iBurlon etaP J201lh the stacking 
of ~200 AGN spectra revealed that the average AGN spectrum is curved in the 15-200 keV 
band. This stacked AGN spectrum can be described using a PEXRAV model with a power- 
law index of 1.8, an energy cut-off of 300 keV and a reflection component due to a medium 
that covers an angle of 2tt at the nuclear source (R~l in the PEXRAV model). These 



- 12 - 



parameters are reported in Tab. 1 of iBurlon et al.l ( 2011'). To c o nvert fluxes from the 15- 
55keV band to the e.g. 14-195 keV band used by iTueller et al.l (|2008[ ) the two factors are 
2.02 and 1.92 (for the power law and PEXRAV model respectively). So we consider the 
uncertainty related to the flux conversion to be of the ~5 % order. 



We cor npare in Tab. P our r esults to those of iTueller et al.l ( l2008l ). ICusumano et al. 



(]2009[ ) and iKrivonos et al.l fl2010h . When comparing INTEGRAL and BAT results one 
has to take into ac count the different no rmalizations of the Crab spectrum that the two 
instruments adopt (IKrivonos et al.ll2010l ). To make a proper comparison we convert the 
INTEGRAL 17-60 keV LogN-LogS to the 15-55 keV BAT band taking into account the 
di fferent normal i zation ^^. It is apparent that there is excellent agreement with the results 
of ITueller et al.l ( 120081 ) both in term o f normalization o f the L ogN- LogS and also in terrn 
of its slope. Both slopes reported by ICusumano et al.l (|2009[ ) and IKrivonos et al.l (l2010l) 



are in agreement with ours, but the density of AGN reported by Cusumano et al. ( 20091 ) is 
sr naller than ours . We fi nd an overall agreement within ~10 % of our results and the ones 
of IKrivonos et al.l (120101 ) as Fig. H] testifies. 



We also c ompa r e the BAT results an d the results of the popula t ion sy nthesis models 
of iGilh et all ( l2007h . iTreister et aP J2OO9I ) and IPraper fc Ballantvnel J2010h . For both the 
Treister et al., a nd Draper &: Ballantyne models we us ed the pred i ctions for the 10-30 keV 
band reported in iBallantyne et al.l (120 111 ) while for the iGilli et al.l ( 20071 ) model we use the 
10-40 keV predictions available onlinqj. Since BAT detect s very few Compton-thick AGN we 



limited the predictions of the models of iGilli et al.l (120071) and 



to objects with LogNn <24. The predictions by ITreister et al. 



Draper fc Ballantvnel (120101 ) 



( 2OO9I ) include objects with 



LogNn >24, however in their modeling the density of Compton-thick AGN is (at BAT 
sensitivities) ~7 % of the total AGN population. The predictions of all models, converted 
to the BAT band using the above prescriptions, are compared to the BAT LogN-LogS in 
Fig. ini It is apparent that there is good agreement with the BAT results at bright fiuxes (i.e. 
> lO^^^erg cm~^ s~^). At the limiting fiux of our analysis (i.e. ~6xl0~^^erg cm~^ s~^) the 
model predictions are beyond the statistical uncertainty of the BAT LogN-LogS as shown 
in Tab. H] and in the inset of Fig. [31 



The model of iGilli et al.l ( 120071 ) is compatible at bright fiuxes with the BAT data, but 



^In order to convert the INTEGRAL data to the BAT band we took into account that 1 BAT-mCrab 
in the 17-60 keV band is 1.22xl0-"erg cm'^ s'^ while 1 INTEGRAL-mCrab is 1.43xl0-"erg cm^^ s-\ 
Subsequently we converted the 17-60 keV fluxes to the 15-55 keV adopting a power law with an index of 2.0. 
This leads to Ff5l'^5=0.878 Flf^ef 



^The lGilli et al.l (|2007() model is available at http://www.bo.astro.it/~gilli/counts.htnil 
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with a steeper slope. On the other hand th ere seems to be a co nsta nt offset between the 
BAT LogN-LogS and t he model predi c tions of iTreister et aLl (120091 ) and lDraper fc Ballantyne 
fl2010h. The model of lTreister et"aD J2OO9I ) reproduces the 14-195 keV BAT LogN-LogS of 



Tueller et al.l ( 120081 ) which is in ve ry good agreement w ith the one publ ished here. However 



a close inspection (see Fig. 1 in ITreister et al.l 120091 ) shows that the iTueller et al.l (120081 ) 
data (reported in ITreister et al.l 120091) have a normalizat ion ~1.6 large r than the original 



measurement reported in iTueller et al.l ( l2008l ). Thus the ITreister et al.l ( l2009l ) is anchored 
to LogN-LogS data with a n ormalization ~ 1 .6 lar ger than observed. This same factor is 
apparent when comparing the ITreister et al.l (|2009[ ) prediction to our data (see Tab. H]). 



The model of iDraper fc Ballantvnel (l2010l) rep r oduce s (see their Fig. 3) the 'correct' 14- 
195 keV LogN-LogS as reported by ITueller et al.l ( 120081 ). However, t heir 10-30 keV mode l 



prediction is at fluxes > 10 ^^erg cm ^ s ^ very similar to the one of lTreister et al.l (|2009[ ) 



and their predicted densities are a factor 1.6-1.8 larger than the measured ones in the 15- 
55 keV band (see Tab. Hj) at the faintest fluxes sampled by our analysis. These findings might 
have some implications for the number of objects predicted to be detected by NuSTAR in 
serendipitous surveys above lOkeV (see §IS])- Given the substantial agreement of the BAT 
and INTEGRAL LogN-LogS (see Fig. Hj), it does not seem likely that the discrepancy in the 
predictions of synthesis models and the >10keV LogN-LogS can be ascribed to a difference 
in the results above lOkeV. 



3.2. The Luminosity Function 



X-ray selected (below lO keV) AGN are known to display a luminosity function that 



Silverman et al. 2008; Aird et al 



evolves with redshift ( see e.g. lUeda et al.l 120031 : iLa Franca et al.ll2005l : iHasinger et al.l 12005 



2010[ ). Our sample reaches a redshift of zf^O.3 where ac- 



cording to t he above wo r ks, the evolution of AGN is significant and can potentially be 
detected. In lAjello et al.l (|2009b( ) we found marginal evidence for the evolution of AGN in 
the local Universe. Here we can make use of our all-sky sample of Seyfert galaxies to test 
this hypothesis. We adopt as a description of the X-ray luminosity function (XLF) a pure 
luminosity evolution model as follows: 



dN 



dVdLx 



^Lxiz),z) = ^Lx/eiz)) 



(3) 



where V is comoving volume element and the evolution is parametrized using the common 
power-law evolutionary factor: 



eu 



:i + z) 



(4) 
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Table 3. Properties of LogN-LogS derived above lOkeV 



Model 


Surface Density^ 


/3 




(10~3deg-2) 





This Work 


6.671°:}! 


2.49in? 


Tueller et al. 2008 


6.5-6.8 


2.42±0.14 


Cusumano et al. 2009 


5.4-5.9 


2.56±0.06 


Krivonos et al. 2010 


7.0-8.1 


2.56±0.10 



^Density at a 15-55 keV flux of 10 erg cm ^ s ^ 



Table 4. Comparison of the BAT LogN-LogS with synthesis models. 



Model 


Surface 


Density^ 




{10-2deg-2) 


This Work 


1.49 


+0.04 
-0.03 


Gilh et al. 2007^ 


2.14- 


-2.35 


Treister et al. 2009^ 


2.31- 


-3.02 


Draper & Ballantyne 2010 


2.52- 


-3.20 



^Density at a 15-55 keV flux of 6x ^^erg 



— 2 —1 

cm s . 



'^To convert the original 10-40 keV counts to 
the BAT 15-55 keV band we have used the fol- 
lowing factors: 0.94 and 1.04 for the power law 
and PEXRAV model. 

•^ To convert the 10-30 ke V counts (reported 



m 



Ballantvne et all 120111 ) to the BAT 15 



55 keV band we have used the following factors: 
1.18 and 1.40 for the power law and PEXRAV 
model. 
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Fig. 3. — LogN-LogS of the BAT AGN (black data points) and best power-law fit in the 15- 
55 keV band. The shaded gray region represents the 1 o uncertainty computed via bootstrap. 
The dashed lines show the predictions of the number counts from the models of Gilli et al. 



(120071 ) , iTreister et al.l ( 120091 ) , and [Draper &: Ballantynd ( l2010l ) . The inset shows a close-up 
view of the distribution at the lowest fluxes. 



We neglect any cut-off in the evolution as this takes place at a redshift that BAT cannot 
constrain (i.e. zs^l see e.g. lAird et al.ll2010l ). For the XLF at redshift zero we use a double 
power-law of the form: 



$(i^x,^ = 0) 



dN 



A 



dLx ln{l(})Lx 



Lx 



71 



+ 



72 



(5) 



A value of k significantly different than zero would point towards an evolution of the 
XLF. In order to the derive the luminosity function of AGN, we adopt the same maximum 
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Fig. 4. — LogN-LogS of the BAT AGN in the 15-55 keV band (filled circles) compared to 
the one derived from INTEGRAL data by iKrivonos et al.l (120101 ) (squares) 



likelihood method described in lAjello et al.l ( l2009bl ). In parti cular we det e rmine the best-fit 
parameters of the XLF by finding the minimum of Eq. 1 1 in lAjello et al.l ( l2009bl ) . 



In order not to include sources which c ould have a non ne gligible contribution to their 
total luminosity from X-ray binaries (see e.g. IVoss fc Aiellolboiol ) and to limit the incomplete- 
ness due to the bias against the detection of the most absorbed sources we derive the XLF 
only for LogLx >41.3. The best-fit parameters for the PLE model in case of evolution [k j^O) 
and in case of no evolution [k =0) are reported in Tab. |5l The first result is that the model 
with no evolution (i.e. Eq. E]) represents a good d escription of th e BAT dataset. The best-fit 
par ameters are c o mpatib le with those reported by lSazonov et al.l (120081 ) . iTueller et al.l (l2008l ) 
and lAiello et~aD J2009bl ). 



Adding the extra k parameter produces only a marginal improvement in the fit. Indeed 
the l og-likelihood improves only by ~5 which corresponds to a ~2(T improvement ( IWilks 
19381 ). This is refiected in the luminosity evolution parameter which is constrained to be 



-17- 



fc=1.38=l=0. 61. The value found here is compatible, but smaller than the value of 2.62±1.18 
reported in lAjello et al.l (l2009bl ) showing that if there is evolution in the XLF of local AGN 
that might be shallower than previously found. We thus believe that the non-evolving XLF 
model is, in view of the marginal improvement of the goodness of fit, a better representation 
of the current dataset. Nevertheless, we will use the PLE model in § |5]to assess the level of 
uncertainty in the prediction for the number of AGN that might be detected by NuSTAR in 
a very near future. 

It is interesting to compare our XLF with that one reported by lSazonov et al. and 



ueller et a 



(120081) in the 17-60 keV and in the 14-195 keV band respectively. 



Sazonov et al. 



+°-^« and 



( 120081 ) and iTueller et al.l (|2008[ ) report a value of the faint-end slope 71 of 0.76_q 20 
0.84^0 22 respectively. The value of our faint-end slope is 0.78±0.08 in agreement with both 
results, but much better constrained because of the larger dataset. When co nverted to our 



band (see p revious section for th e conversion factors), the break luminosity of ISazonov et al. 



(2008 


) and 


Tueller et al. 


(2008) 



h2.0 



X 10 



43 



erg s" 



-1 and 3.7tls x 10^^ erg 



s ^ while we measure 5.1±1.4 x 10^^ erg s ^ again compatible, but better constrained. 



In orde r to display t 



2005 



(11997) and Mivaii et al. 



leLF we rely on the "No^«/N'»'^'" method devised by 
( 200l[ ) and employed in several recent works (e.g. 



La Franca & Cristiani 


La Franca et al. 





Hasinger et al.ll2005[ ). Once a best-fit function for the LF has been found, it is possible 



to determine the value of the observed LF in a given bin of luminosity and redshift: 

j^obs 



X,i, ■ 



mdl 



(6) 



where Lx,i and Zi are the luminosity and redshift of the i*'^ bin, ^"^'^''{Lx,i, Zi) is the best-fit 
LF model and N"^'^ and AJ"*^^ are the observed and the predicted number of AGN in that 
bin. 

Fig. [5] shows the best- fit non-evolving mod el (i.e. fc=0) in comparison with the XLF of 
Sazonov et al.l (120081 ) and iTueller et al.l (120081 ). In general there is very good agreement 
between the XLFs derived in all these works. Our results are als o in agreement wit h 
those obtained in the <10 keV band. I ndeed , for the bright end slope iBarger et al.l (120051 ). 



La Franca et al.l (120051 ) and lAird et all (I2OIOI ) obtain 2.2± 0.5, 2.36 



- tnl? and 2.55±0. 12 respec- 



tively while we measure 2.39± 0.12. For the faint-end slope iBarger et al.l ( l2005l ) , iLa Franca et al. 
J2005I ) and lAird et aP J2010[ ) report 0.42±0.06, 0.97lJ]:[J? and 0.58±0.04B while we measure 



7i =0.78±0.08. The agreement for the faint-end slope is not as a good as for the bright-end 
slope and there is some scatter (that appears to be systematic in origin) in the <10keV 



' Aird et al.l |2O10l ) reports also 0.70±0.03 for their PLE model. 
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measurement while there is substantial agreement above lOkeV (see above discussion). An 
excessive flatness of the faint-end slope might be linked to the role of Compton-thick AGN 
which are difficult to detect because the ab sorption push e s thei r observed luminosity below 
the survey threshold. Indeed, as shown in iBurlon et al.l ( 120111 ) 2-10 keV surveys are more 
biased in the detection of Compton-thick AGN than surveys above lOkeV. Depending on 
the redshift distribution of the sources (and thus on the k-correction) the bias might not be 
the same for different 2-10 keV surveys. 

For completeness we also report in Tab. the best-flt parameters to the XLF of all 
AGN excluding the Compton-thick ones of TabJH Because Compton-thick objects represent 
a small fraction of the AGN detected by BAT, there is very little difference between the 
XLF of all AGN and that of AGN with LogN^ <24. However, in the next section this XLF 
will be useful to compare the space density of Compton-thick AGN to that of Compton-thin 
(LogNj^ <24) AGN. 



3.3. The Space Density of Compton-thick AGN 



We derive the space density of Compton-th ick AGN with 24<LogNH <25 using the 
1/Vmax non-parametric method ( jSchmidtl Il968l ). We de-absorb the luminosities of the 
BAT Compton-thic k sources reported in Tab. [2] using the correction function reported in 
Bur Ion et al.l ( 1201 ll ) (see their Fig. 11). This correction function was derived for the average 
properties of the Compton-thick AGN detecte d by BAT and takes into a ccount photoelectric 
absorption as well as Compton scattering (see iMurphy fc Yaqoobll2009l . for details). 



Fig. [6] shows the luminosity function of the Compton-thick AGN detected by BAT while 
Fig.[7]reports the cumulative space density. In this latter case the uncertainty were computed 
via bootstrap with replacement. We find that the space density of Compton-thick AGN with 
a de-absorbed luminosity greater than 2 x 10*^^ erg s 



Mpc- 

absorbed luminosity of 10 erg s~ the de n sity becomes 2.1+^;^ x 10"^ Mpc" 



is r.dtti X 10"^ 

-1.6 



Above a de- 
As shown in 

Fig. [H] the model predictions of iGilli et al.l (120071) are com patib le within the statistic a l error 
with our space density estimates. iTreister et al.l (120091 ) and iDraper fc Ballantynd (120101 ) 
estimated a space density of LogLx >43 erg s~^ sources of respectively 2.2tli x lO'^Mpc"^ 
and 3-7x10^^ Mpc~'^. These estimate are below ours, but compatible within 2a and la 
respectively. 

In the same Fig. [6] we plot also the space density of AGN with LogNn <24 (from the 
previous section). While there is substantial agreement between the two, our analysis seems 
to point to the fact that the space density of Compton-thick AGN is larger than the one of 
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Fig. 5. — Luminosity function in the 1 5-55 keV ban d of th e BA T AGN (non evolving: model) 
compared with the measurements of iTueller et al.l ( l2008l ) and ISazonov et al.l ( l2008l ). The 
data point at a luminosity < 2 x 10^^ erg s~^ was not fitted to avoid problems related to 
incompleteness in detecting Compton-thick AGN and to avoid cor itamination from so urces 
whose flux might be dominated by the emission of X-ray binaries (IVoss fc Ajelldl2010l ). 



all other classes of AGN. By allowing the normalization A of the XLF of LogNn <24 AGN 
to vary we find that ALogNH>24=l-4xy4LogNH<24: i-e. the space density of Compton-thick 
AGN is 1.4 times larger than that of LogNn <24 AGN. A similar results was obtained by 
Delia Ceca et al.l ( l2008bl ) who derived indirectly the space density of Compton-thick AGN 
as a difference between the density of optically selected AGN and that of X-ray selected 
AGN with LogNn <24 (see the aforementioned paper for more details). In their study 
they also find that the space density of Compton-thick AGN with LogLx >43erg s~^ is 
~1.6x 10^^ Mpc^'^ which is in good agreement with the value found here. 



10-^- 



10- 



10" 



Space density of Compton-thick AGN 



Predicted space density (Gilli et al. 2007) 



Space density of Compton-thin AGN 

J I I I I I I I I I I I I I I I I IM 1 




10 



10-^ 

Lx [10""* erg s'^] 



10 



Fig. 6. — Spa ce densities of C ompton-thick AGN detected by BAT compared to model 
prediction from lGilli et al.l ( 120071 ) and the space density of Compton-thin (LogN n < 24) AGN 
from ^ 13.21 Luminosities were de-absorbed following the method outlined in iBurlon et al. 



torn 



4. The Anisotropic Local Universe 

The BAT sample represents an incredible resource to study the properties of the nearby 
Universe as it is a truly local sample of AGN. Indeed, the median redshift of the radio- 
quiet AGN detected by BAT is ~0.03 which corresponds to ~120Mpc. On these scales the 



structure of the local Universe is known to be in 



the matter density distribution (ICappelluti et al. 



lomogeneous and AGN are known to trace 
2010). The BAT AGN can thus be used to 



trace the large-local scale structure. Krivonos et al. " (j2007 ) reported, using 68 AGN detected 



by INTEGRAL, an anisotropy of the spatial distribution of local AGN. With an AGN sample 
~6 times larger than the INTEGRAL sample it is possible to identify the anisotropic spatial 
distribution of AGN to a higher precision than previously obtained. 
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Fig. 7. — Cumulative space densities of Compton-thick AGN detected by BAT (thick red 
solid line) with 1 a confidence contours (dashed line) generated from the analysis of the 
bootstrapped samples (thin gray lines). 

The following approach is adopted in order to assess the anisotropy of the spatial distri- 
bution of AGN. We first separate the sample in two redshift bins z<0.02 and 0.02<z<0.15 
corresponding to a distance smaller (or grater) than 85 Mpc. The choice is dictated by the 
fact that the largest contrast in AGN density is expected to be observed nearby from us. 

For each of these sub-samples we create a LogN-LogS that yields the average surface 
density of AGN for that redshift slice. Taking into account how the BAT sensitivity varies 
across the sky we generate (for each redshift bin) 1000 realizations of AGN sets which are 
isotropically distributed in the sky. Then for each direction in the sky we count how many 
AGN BAT has detected within a radius of 20° (typically this number oscillates around 
10) versus the expected number for the isotropic case. For each of these sky positions we 
compute the fractional over-density of AGN: i.e. the ratio between the number of detected 
AGN and the number of expected objects if AGN were isotropically distributed. From the 
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1000 realizations we also compute the error (and the significance) connected to the fractional 
over- density. 

Fig. [S] and Fig. 12] show how the fractional over-density of AGN (i.e. the ratio of detected 
AGN to the expected number for the isotropic case) and its significance change across the 
sky for the lower redshift bin. It is apparent that there is a marked contrast in the AGN 
density within 85Mpc, with the AGN density varying by a factor of >10 across the sky. 
The most significant over-density is seen in the first redshift slice (i.e. within 85 Mpc) with 
prominent structures significant at >4cr. 

There is a clear over-density of AGN (a factor >4 larger than the isotropic expectation) 
in the direction of the super-galactic plane. In particular the most dense region is found 
to be at / = —54° extending from b =~ 15° to b =~ 50°. This corresponds to the posi- 
tion of the Hydra- Centaur us super-cluster (z=0.01-0.02) which represents one of the largest 
structures in the local Universe. The second most dense region (connecting to the south 
of the Hydra- Cent aurus super-cluster) can likely be identified with the 'Great Attractor' 



( iLynden-Bell et al.lllQSSi ). In Fig. [H]the positions of the Great Attractor, the Centaurus and 
Hydra superclusters are shown with squares (from bottom left to upper right, respectively). 
The redshifts of the BAT AGN in these dense regions are in agreement with the redshifts 
of these massive s tructu res. Our results appear to be in good agreement with the ones 



of iKrivonos et al.l (120071 ). but the improved statistics allow us to locate more precisely the 



over-density of AGN in the nearby Universe. 

In the 0.02<z<0.15 the anisotropy is less pronounced and less significant as well, with 
the most significant structures being < 3 cr. A comparison with our simulations (which 
allow us to account for the trial factor) shows that the 0.02<z<0.15 over-density map is 
indistinguishable (with the present data) from the isotropic expectation. How 'ever, the loca l 



environment is known to be highly inhomogeneous up to (at least) zf»0.05 ( jJarrettl 120041 ). 
The reason why BAT does not detect this anisotropy in the local scale structure is ultimately 
due to statistics. Indeed, the typical clustering length of the BAT AGN is known to be 5- 



8 Mpc ( Cappelluti et al. 2010 ). Within 85 Mpc, our 20 degree search cone subtend^ a length 



of ~16 Mpc which is comparable to the clustering length of AGN. However, for the average 
redshift (~0.05) of the 0.02<z<0.15 sample, the 20 degrees search cone subtends a length of 
70 Mpc which is ten times larger than the typical clustering length of AGN. This contributes 
to dilute any over-density signal. In order to have a similar resolution as for the <85 Mpc 
sample, one should adopt a search cone of 5 degrees. However, with the current statistics we 
would expect less than 1 AGN in such cone. 



'''The average redshift of the z<0.02 sample is 0.01. 
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Finally we checked if the anisotropy depends on any of the AGN parameters (i.e. flux, 
luminosity, type, etc.). To this extent we isolated the AGN in the direction of the Hydra- 
Centaurus super-clusters and of the Great Attractor and created a LogN-LogS and a lumi- 
nosity function. The first result is that the XLF of these AGN is in good agreement with 
that of the entire population reported in Fig. |5l The LogN-LogS of the AGN in the direction 
of the super-clusters is shown, in comparison with the LogN-LogS of the whole population, 
in Fig. [TOl It is clear that there is an excess at bright fluxes. This finding in conjunction 
with the fact that the XLF does not change indicates that the sources that contribute to 
the over-density are the brightest sources (i.e. the proximity of these local structures makes 
these sources appear with a bright flux). Indeed, the median redshift of the AGN in the di- 
rection of the super-clusters (z^O.015) is markedly smaller than that of the whole population 
(zRiO.03). 

It is also interesting to note that the fraction of Sy2 galaxies in the direction of the 
super-clusters is larger than average. Indeed, the fract i on of objects classified as Sy2 in 
our whole sample and in the sample of ICusumano et al.l (120101 ) is ~34 % of the total AGN 
population. This ratio would increase to ~45% if all the objects with no opictal classification 



reported in this or the ICusumano et al.l (120101 ) sample would be in reality Sy2s. While it is 
reasonable to expect that more than 50 % of the objects lacking an optical classification are 
Sy2 galaxies, it is unlikely that all of them are Sy2s. So this represents an extreme scenario. 



The fraction of Sy2 object becomes ^ 
the super-clusters and increases to 58 ±8 ' 



■^50 % if we consider only the AGN in the direction of 
^ if we restrict to sources within 85 Mpc. This seems 



in agreement with what reported by iPetrosiaru (119821 ) that Syls are more often found to be 
isolated than Sy2s. In dense environments encounters between galaxies produce gravitational 
interactions that can trigger gas infiow towards the central hole and produce AGN activity. 
Galaxy merging is known t o provide an efficient way to funnel large amount of gas and dust to 



the c entral black hole (e.g. iKauffmann fc Haehneltll2000l : IWyithe fc Loebll2003l : ICroton et al. 



2006h . Recently, IKoss et al.l J201oh found that 24% of all the hosts of the BAT selected 
AGN have a companion galaxy within 30kpc. This suggests that for a fraction of moderate 
luminosity AGN merging is a viable triggering mechanisms. However, in dense environments 
major merging is not the only process that might be at work. Indeed, galaxy harassments 
(i.e. high spe ed encounters be tween galaxies) can also drive most of the galaxy's gas to the 
inner 500 pc (ILake et al.lll998l ) and thus trigger AGN activity. Also minor merging, where 
the ratio of the masses of the merging galaxies is >3 ( and probably around ~10), can lead to 
Seyfert level of accretion (IHopkins fc Hernquistll2009l ). With large quantities of gas available 
around the hole, in the early phase of accretion the AGN might have a higher probability of 
being identified as a type-2 AGN. This might explain the over-abundance of Sy2 galaxies in 
dense environments. 
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Galactic Longitude 



Fig. 8. — Fractional over-density of AGN with respect to the isotropic expectation within 
85 Mpc. The color scale shows the ratio of BAT-detected AGN (within 20 degree cones) to the 
average number of AGN expected in the same area if sources were isotropically distributed. 
The map shows the regions where this ratio is the largest (red) and where it is the smallest 
(blue). The squares show the approximate position of the most prominent super-clusters 
(see text for details). 



5. Predictions for NuSTAR 



In this section we provide predictions for the number of objects that NuSTAR might 
see in different types of blank-field surveys. The predicted number counts are obtained by 
extrapolating the BAT LogN-LogS of § 13.11 to lower fluxes under the assumption that the 
slope of the LogN-LogS does not change. Indeed it is reasonable to expect so since there are 
strong indications both from observations in the 2-10 keV band and from modeling that the 
source count distribution breaks only at fluxes < 10^^^ erg cm~^ s~^. 

The source count distribution of AGN is ve ry well determined in the 2-10 keV ban 
down to fluxes of < 10~^^erg cm~^ s~^ (see e.g. iRosati et al.l 120021 : ICappelluti et al.l 120071 : 



^For a typical AGN spectrum the 2-10 keV flux is on average 20% larger than the 15-55 keV flux. 
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Galactic Longitude 



Fig. 9. — Significance of tlie density features of Fig. |H] expressed in number of a. 



Xue et al. 


2011 


). Its slope is 


For example 


Cappelluti et al. 



10 



-14 —2 - 

erg cm s 



results. A similar result is found by iRosati et al. 
predict a break in the LogN-Log S at < 10 
2OO9I : IPraper fc Ballantvii3l201oh . 



fl2002f ) 



1 4 —2 

erg cm s 



Al so population sy i ithesis models 
-1 Jcilli et al.l[2007l : IXreister et al. 



The extrapolation to fluxes a factor >50 fainter than those sample by BAT necessarily 
introduces some uncertainty related to neglecting the evolution of AGN. This uncertainty 
will be gauged later on in this section. However, the fact that BAT and NuSTAR sample 
(almost) the same energy band removes other sources of uncertainties. NuSTAR will likely 
perform three different types of surveys: a shallow, a medium and a deep survey. The shallow 
survey, performed combining short (5-10 ks) exposures, might extend over ~3deg^ reaching 
fiuxeqj ~ 10~^^erg cm~^ s~^. The medium survey will likely cover ~ldeg^ with ~50ks 
pointings reaching a 10-30 keV flux of ~ 5 x lO^^'^erg cm^^ s^^ while the deep survey is 



^Typical fluxes for NuSTAR are quoted for the 10-30 keV band. Here we have converted those fluxes to 
the 15-55 keV band adopting a power law with a photon index of 2.0. 
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Fig. 10. — LogN-LogS (15-55 keV) of the AGN in the direction of the Hydra-Centaurus 
super-clusters and of the Great Attr actor compared to the LogN-LogS of all AGN. 



expected to reach a 10-30 keV flux of 
an area of ~0.3deg^. 



2 X 10 



-14 —2 

erg cm s 



(using 200 ks pointings) over 



Ballantyne et al.l (120 111 ) reported the number of expected AGN, in NuSTAR surveys, 



as predicted from different population synthesis models. Since the NuSTAR field of view 
is smaller than the surveyed area, NuSTAR will have to use a tiling stra tegy. Two tiling 



strat egies can be foreseen: a corner shift and half-shift survey (see also iBallantyne et al. 



20111 ). In the corner-shift strategy, the survey area is covered by non-overlapping pointing. In 
the half-shift strategy the distance between pointings is half the size of the FOV. The corner- 
shift strategy reache s deeper fluxes w h ile th e half-shift reaches a more uniform exposure of 
the surveyed area. IBallantyne et al.l ( 120111 ) concluded that the half-shift strategy yields 
the larger number of AGN. Thus, in or d er to make a proper comparison we adopt the sky 
coverages reported by IBallantyne et al.l (120111 ) for the half-shift tiling strategy for the 10- 
SOkeV band. Converting the sky coverages from the 10-30 keV to the 15-55 keV is rather 
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easy and almost error free. Indeed for the power-law and the PEXRAV model discussed 
in § 13.11 we get a conversion factor of 1.16 and 1.18 respectively. Thus using a LogN-LogS 
derived in an overlapping band reduces the uncertainties due to flux conversion to less than 



The number of objects predicted, extrapolating the LogN-LogS, for 4 different survey 
fields are reported i n Ta b. [6] and they are compared to the pre dictions of the models o f 
Treister et al.l (120091 ) and iDraper fc Ballantynd (120101 ) reported in iBallantyne et al.l ( 120111 ) . 
It is clear that our predictions lie substantially lower than the ones of population synthesis 
models. The left panel of Fig. [11] shows that (at NuSTAR sensitivities) the predictions of 
synthesis models lie a factor ~3 above the extrapolation of the LogN-LogS of BAT AGN. 

This disagreement is likely due to two reasons. Part of it is due to the fact that synthe- 
sis models lie systematically above the BAT LogN-LogS as shown in § 13. 1[ The second one 
is that the simple extrapolation of the BAT Log N-LogS is not able to capture the evolu- 



tion of AGN which i s rather well established (e.g. iMiyaji et al.l 1200 ll : iLa Franca et al.l 12005 
Hasinger et al.ll2005l ) and expected also in the >10keV band. 



In order to correct for the first problem ( i.e. the over-pred i ction of AGN i n the BAT 
band) we arbitrarily renorn ialize the models of iGilli et al.l (|2007| ). iTreister et al.l (120091 ) . and 
Draper &: Ballantynd (120101 ) to fit the BAT data at bright fluxes and then we convolve them 
with the same sky coverages described above. The predictions for the 're-normalized' models 
are reported in the lower part of Tab. [H] and are on average compatible within ~1 a with the 
extrapolation of the BAT LogN-LogS (as also visible in the right panel of Fig. [TTj) . This is 
reflected into a lower number of predicted AGN detections for NuSTAR in the 10-30 keV 
band. 

In order to correct (or to gauge the uncertainty) due to neglecting any evolutionary 
effect in the LogN-LogS we rely on the best-fit PLE model of § 13. 2[ Even if marginally 
significant we take the best-fit parameters at their face value and derive the number of 
expected objects in NuSTAR fields. These are reported in the last row of Tab. [61 As it 
can be seen the prediction from the evolutionary XLF are in fairly good agreement with the 
predictions from the 'normalized' synthesis models. We thus believe this set of predictions 
(i.e. lower part of Tab. [6]) for the number of AGN detectable is realistic. However, the 
ultimate number of detected AGN will likely depend on how the sensitivity will vary across 
the survey fields and, for the smallest fields, on cosmic variance. 
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Table 5. Parameters of fitted Luminosity Functions. Parameters without an error 
estimate were kept fixed during the fitting stage. 

Model 7i 72 L^^ k 

No evolution 113.1±6.0^ 0.79±0.08 2.39±0.12 0.51±0.14 
PLE 109.6±5.3^ 0.78±0.07 2.60±0.20 0.49±0.10 1.38±0.61 

Noevol., noCT 122.4±6.6^ 0.72±0.09 2.37±0.12 0.48±0.13 

''In unit of lO^'^Mpc"^. 
'^In unit of 10^^ erg s~^. 



Table 6. Predictions f or NuSTAR surveys adopting the half-shift coverages reported in 



Ballantyne et al.l ( 1201 ll ). The lower part of the table shows the prediction of the models 



renormalized to match the BAT LogN-LogS. 



Model 


BOOTES 


COSMOS 


ECDFS 


GOODS 




(9.3deg2) 


(2dcg2) 


(0.25 deg2) 


(0.089 deg2) 


This work (no evo.) 


47+16 
^' -13 


01+13 
'^-'--10 


20t^° 


16t^ 


Draper & Ballantyne 2010 


126 


91 


62 


51 


Treister et al. 2009 


107 


77 


52 


42 


Draper & Ballantyne 2010 


66 


46 


31 


26 


Treister et al. 2009 


68 


51 


33 


27 


GiUi et al. 2007 


94 


60 


39 


25 


This work (evo.) 


61 


42 


28 


23 
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Fig. 11. — LogN-LogS of the BAT AGN (black data points) and best power-law fit ex- 
trapolated to the sensitivity that NuSTAR will achieve. The left panel shows the original 
prediction of synthesis models, while in the right panel the predictions have been renormal- 
ized to match the AGN densities measured by BAT. In both panels the vertical dashed line 
shows the sensitivity reached by NuSTAR in a short pointing (i.e. 5-10 ks). 

6. Summary and Conclusions 

The analysis presented here shows the power of an all-sky survey at hard X-rays to study 
the AGN in our local environment. Thanks to the large field of view and high sensitivity 
BAT has detected 428 AGN all-sky above lOkeV with negligible (<5%) incompleteness. 
This represents the largest complete sample of AGN detected so far. Below we summarize 
our findings. 

• The BAT AGN sample spans 10 decades in luminosity comprising objects detected at 
distance of ~lMpc up to redshift ~3.5. The AGN sample can be divided into objects 
whose emission is dominate by the accretion disk/corona (i.e. Seyfert galaxies) and 
jet-dominated sources (i.e. blazars and radio galaxies). Seyfert galaxies are detected 
at low redshifts and low luminosities while jet-dominated sources, that are ~15 % of 
the whole sample, are detected at high luminosity and high redshift. 

• Samples of AGN detected above lOkeV are instrumental to determine the size of 
the population of Compton-thick AGN that are still detected in very little numbers. 
While a detailed measurement of the absorbing column density of all the AGN in the 
BAT sample is left to a future publication, we cross-correlated our sample with known 



- 30 - 



catalogs of bona-fide C ompton-thick AGN. The BAT sample already comprises 15 (out 
of the 18 reported by iDella Ceca et al.l l2008al ) bona-fide Compton-thick AGN and 3 
likely candidates. The observeci^ fraction of Compton-thick AGN, relative to the whole 
population, is thus ~5%. We showed that BAT will likely not detect the rest of the 
known candidate Compton-thick AGN. Since the BAT sensitivity still improves with 
time, future AGN samples detected by BAT will likely contain previously unstudied 
Compton-thick AGN. A few (1 or 2 objects) might be present already in this sample. 

We performed a robust analysis, using bootstrapping, of the source counts distribution 
of the Seyfert-like objects detected by BAT. The BAT LogN-LogS is consistent with 
Euclidean down to the lowest fluxes spanned by this analysis (i .e. 6xl0~^^ erg cm~ ^ 



The agreement between our and the INTEGRAL results (iKrivonos et al.l 120101 ) 



shows that the LogN-LogS of AGN selected above 10 keV is establis hed to a precision 
of 10%. The population synthesis models tha t we have tested (i.e. iGilli et al.l 120071 : 
Treister et al.ll2009t iDraper fc Ballantynell2010[ ) are able to reproduce the BAT LogN- 
LogS at the brightest fluxes, but overestimate it at the lowest fluxes spanned by this 
analysis. 

We derived the luminosity function of (local) AGN and tested for its possible evolu- 
tion with redshift. Even with our large sample of AGN the evidence for the evolu- 
tion of the XLF are at best marginal (i.e. ~2cr). The BAT data are well described 
by a non-evolving XLF which is modeled as a standard double power law. We find 
that the slope of the faint end is 71 =0.79±0.08 while that on e of the bright end is 
791= 2. 39±0. 12. These va lues are in good agreement with those of ISazonov et al.l (l2008[ l 
and lTueller et al.l ( l2008l ). but better constrained. Given the small FOV, NuSTAR will 
not be able to directly constrain the properties of the low-luminosities low-redshift 
AGN. The BAT sample and the BAT XLF are thus instrumental to determine, in 
connection with the future NuSTAR samples, the evolution and growth of AGN. 

We derived the luminosity function of Compton-thick AGN and found that at redshift 
zero their space density is 7.9l2'9 ^ 10~^ Mpc~^ for objects with a de-absorbed lumi- 
nosity larger than 2x10^^ erg s~^. Our measurement is slightly larger, but compatible 
within uncertainties, with the prediction of synthesis models. 

The BAT samples of Seyfert galaxies is truly a local sample (median redshift 0.03) 
and can be used to study the spatial distribution of AGN in the local Universe. We 



See lBurlon et all ([201 ll ) for the bias that also instruments above 10 keV have in detecting Compton-thick 
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detected significant over-density features in the spatial distribution of AGN located 
within 85 Mpc. The densest regions show a density of AGN that is up to ~5 times 
larger than the average all-sky density. These dense regions can be identified with 
the most prominent nearby super-clusters: i.e. the Hydra-Centaurus super-cluster and 
the 'Great Attractor'. The fraction of Sy2 galaxies (with respect to the total AGN 
population) appear to be larger than average in the direction of the over-dense regions. 
This evidence might support an evolutionary link where close encounters of galaxies 
trigger AGN activity whose first appearance is obscured by dust and gas. 

• The BAT and the NuSTAR energy bands overlap and it is thus possible to derive 
straightforward predictions, from the LogN-LogS or the XLF, for the number of AGN 
that NuSTAR might detect in survey fields in the near future. We find substantial 
agreement in the number of predicted objects if the predictions of population synthesis 
models are renormalized to match, at the lowest fluxes, the BAT LogN-LogS. 

Owing to the capability of detecting the local (those within ~200 Mpc) , relatively low- 
luminosity AGN, but also high-redshift blazars, the all-sky hard X-ray surveys - such as the 
Swift-BAT survey discussed in this paper - have the unique potential to study simultaneously 
both the nearby and the early Universe. Such hard X-ray surveys represent a unique resource, 
since for many years, they will remain the most sensitive probes of the accretion history in 
the Universe. 

MA acknowledges extensive discussions with Nico Cappelluti. MA acknowledges Eze- 
quiel Treister and David Ballantyne for making the predictions from their models available 
in a machine readable format. The authors acknowledge the comments of the referee which 
helped improving this paper. 
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